Self-rotating wave approximation via symmetric ordering of ladder operators 
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We show how some Hamiltonians may be approximated using rotating wave approximation meth- 
ods. In order to achieve this we use the algebra of boson ladder operators, and transformation for- 
mulas between normal and symmetric ordering of the operators. The method presented is studied 
in two special cases; the Morse and the Mathiue models. The connection with regular perturbation 
theory is given and the validity of the approximation is discussed. 

I. INTRODUCTION 



In this contribution we produce approximations for some Hamiltonians such as the ones related with the Morse 
potential and cosine potential (quantum rotor, or Mathiue differential equation). The approximation is based on the 
rotating wave approximation (RWA) [l|, 0| and symmetric ordering of ladder operators "B . The RWA is widely used 
in quantum optical systems and particularly in Jaynes-Cummings types of models [U, H[ and its extensions jy] , for 
example in ion traps 0, cavity quantum electrodynamics 0] and population transfer in atoms and molecules 0]. 
Also the breakdown of the approximations in similar systems has been considered in several papers [lfj. The RWA 
is usually used for two (or more) interacting subsystems, and in all the above examples, interaction between different 
subsystems are considered. When passing to the interaction picture with respect to the free Hamiltonian, one remains 
with a time dependent Hamiltonian that under certain conditions on the frequencies involved may take a simpler 
form by neglecting terms that oscillate rapidly @, 0|- A phenomelogical motivation for the approximation is that 
the neglected terms usually describe non-energy conserving processes, which in the Jaynes-Cummings model discribe 
simultaneous excitation of the atom/ion and the quantized field. The neglected terms, usually very small within the 
axperimental parameter regimes, give rise to the Bloch-Siegert shift of the energies In this paper we show that 
we can also apply a RWA in a single Hamiltonian systems. Developing the potential in a Taylor series and grouping 
the quadratic term of the Hamiltonian with the kinetic energy to produce a harmonic oscillator term plus an infinite 
sum. The frequency of the artificially produced harmonic oscillator (HO) is then used as a reference to apply the 
RWA to the remaining sum. We call this HO a self-HO for the Hamiltonian studied. Using the underlying algebra of 
the oscillator ladder operators and the RWA, closed forms of the infintie sums are given in two special cases. Only 
Hamiltonain systems containing bound or quasi bound states may be approximated, and the validity depends on the 
relative "deepth/ width" of the potential; the deeper and narrower potential the better approximation. Thus, it is 
best suited for low excited bound states, which is verified by numerical simulations. It is also shown that the RWA 
gives the result of regular first order perturbation theory of the self-HO, which again underlines the validity regimes 
of the approximation. 

We proceed as follows, in the next Section we introduce the self-HO and its Fock states of an arbritary Hamiltonian 
and derive expressions for the diagonal elements of this Hamiltonian in the Fock basis. This is achieved using the 
properties of the ladder operators. In sections ITTTl and lTVl the connection with the RWA is given and two specific models 
are considered; the Mathieu and Morse equations respectively. The validity of the approximation is investigated in 
Section|V]for both studied cases. Secton lVll is left for a summery, and in the appendix lAl is is shown how the methods 
may be used in a more genaral way for calculating various sums. 



II. DIAGONAL MATRIX ELEMENTS IN THE FOCK BASIS OF AN ARBITRARY HAMILTONIAN 

Given an arbitrary Hamiltonian we expand the potential term in a Taylor series (in the following we consider for 
simplicity the case of unit mass and h = 1) 

H = £ + V(x) = V(»)(0) + VW(0)i + £ + + £ (1) 

k=3 

with uj 2 — V^ 2 )(Q). We let this frequency to define the ladder operators [l2j for the Hamiltonian systems as 
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The eigenstates of the number operator N = a^a is the Fock states, labeled \N). The ladder operators incrase or 
decrease the quantum number N when acted on a Fock state; a^N) — \/ N + T|iV+ 1) and a\N) = \/N\n— 1). Thus, 
only terms of the infinite sum 

5 = E^74r( at + a ) • ( 3 ) 

fc=3 

of equation ([T]) that contain an equal number of aJis and a:s will contribute to the diagonal matrix elements Sn — 
(N\S\N). When considering those elemenst we may therefor neglect all other terms of the sum. We have 

I ( h/o J : (a t a) fc/2 '-w, k even 

(at+a) fc ^«j v fc / 2 ; (4) 

k odd 

where : (a^a) fe :w denotes the Weyl (symmetric) ordering of the operator N k . It may be transformed into normal 
ordering using [131 ] 

■■^af: w ^J2ll( k l ) 2 ^ k - l - k ~ l - (5) 
i=o ^ ' 

Thus, the new sum S(x,p), when the "off-diagonal" terms have been neglected (this sum will now depend on the 
operator p), may be written as 



2 fe (fc!) 2 2 l V I . 

fc=0 Z=0 v ; v 7 

The second sum may be taken to infinity (as we may only add zeros) 



1=0 k=0 y ' 



For k < I the above expression is zero and thus we make the substitution j = k — I and get 

1 ~ V (2j+2l) 

2^V.^ 2i(j)\ 



- " 1 A V^)(0) ^ 

i=0 j=0 



By using the identity 



+ . . iV! 

abaJ = _^ ( 9 ) 

the sum is diagonalized in the number basis {|iV)} and we get the diagonal elements of (JTJ) 

{ N\H\N) = „(N + -) + g ^ g 2J( ^ (10) 

where we have used that (iV|^ 2 /2|iV) = (A|w 2 x 2 /2|A) = |(iV + |). Note that the Fock basis used to calculate the 
diagonal elements of the Hamiltonian are obtained from the quadratic term of the Hamiltonian and we call them \N) 
instead of the usual notation \n). In the following sections we apply the above to some specific examples. 



III. COSINE POTENTIAL 



Using the fact that we can obtain closed forms for the diagonal elements of some of the potentials studied in Section 
im and noting that the RWA is used to keep constant terms (terms that do not rotate), in this Section we show that 
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for some Hamiltonians we can produce approximations under certain circumstances, namely, when some parameters 
allow us to perform the RWA. 

One very important equation for scientists is the Mathieu equation [3, EH, EH] which is identical to the Schrodinger 
equation with a sin or cosin potential, also known as the quantum rotor. Durning the last decade it has gained a 
new shove of attention due to the growing field of cold atoms in optical lattice s 11711. App roximate results of the 
eigenf unctions and eigenvalues of the Mathieu equation have been studied earlier [lSL 1 19L l2fj • Known approximation 
methods from physics, such as the WKB and the Raman-Nath have been applied to the Mathieu equation [li|. Given 
the Hamiltonian (reminding that we set h = 1 and m = 1) 



P 2 

H=—-glcosqx (11) 
we can expand the cosine as (we neglect the constat term as it only displaces the energies) 

H= P - + ^-glf { -p^(-l) k . (12) 
2 2 Z^i ( 2fc ) ! 

With the above definition of the creation and annihilation operators © , we get 

[go<l~ . . P -+ fdoQ- ■ P rt9\ 

a = \j- x + l VW q > 4t = fT 8 -^' (13) 



and rewrite the Hamiltonian as 



fc=2 



(14) 



In the case go 3> | we can do RWA, and thus go to a rotating frame with respect to the free Hamiltonian, U(t) = 
exp (— igoqaJat). Using the relations 

U\t)aU{t) = ae i909t , rf(t)tfU(t) = ate"* 909 ', (15) 

and by only keeping non-rotating tyerms in the Hamiltonian we find 

~ (-1) 



H = 9oq (N + l)- gl g {£) t± : (atfi)* : W (16) 
We now proced as in section [III Inserting ([5|) into (fT6|) we have H = goq (^N — 1/2^ — g^S with 



,_,,_„ ,W (*!) 2 2'V», 
Rearranging terms and like in equation |7]) we can sum I to infinity to get 

oo oo 



By expressing 



we can write S = Si + 52 with 

OO OO / \ ft 



°° 1 °° / a \ ft f_ 



2 l l\f^\2g J (fc-Z)! 2 
=o fc=0 v yu/ v ; 



(19) 



(20) 



4 



and 



&= (4 t a (21) 

Note that, like previous in section [III the second sum in equation ([20]) may be started at k = I, i.e. 



DO ^ DO 



or 



2 l ll^\2g J {k~l)\ 2 ' 



n 1 / 



or using that a^^a? — 



Nl 



(N-j)\ 



Si = : J ( ^-tfa ) : (23) 



S x - e~^L^ ( A- ) (24) 
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The result obtained here agree with previous results were approximations to the Mathieu equation are obtained [l£ 
The result also gives an interesting relation between zeroth Bessel function and Laguerre polynomials. 
Combining the results we find the diagonalized RWA Hamiltonian 

9oq ( ' tCt , 1 



s - ! n' r+ 5r*" ,B M£j- (25) 

If we develop the Laguerre polynomials in powers of q/ga, which should be valid within the regimes of the RWA, and 
remain to second order we obtain 

H^ gm {N+ l -)-^(N 2 + N+ l -)-gl (26) 

It should also be pointed out that the method also works for sup erimposed cosine potentials; qj cos(tfex). This 
has interesting applications in for example solid state physics [2l| and cold atoms in optical superlattices [22j . 

We conclude this section by making an analogy betwen the RWA and first order perturbation theory. Given the 
Hamiltonian 

H = g q (n + (27) 

with eigenstates \N), we perturb it with 

V = -g 2 cos(qx)-^fx 2 . (28) 

The first order corerection to the energy [23[ becomes 

5E = (¥„M* B ) = -g 2 e-*kL n - if (n + ~) , (29) 

which regains exactly the same result as the one obtained from the RWA method. This relation between first order 
perturbation theory and the RWA method is easily shown, using Eq. (|A4j) . to be valid in any general case. 
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IV. MORSE POTENTIAL 

Now let us look at the Morse potential 0, HH , from which we get the Hamiltonian 

fi 2 

H = I -j + \ 2 {l-cxp[-a(x-b)}) 2 (30) 

This kind of Hamiltonian is commonly used to describe properties of diatomic molecules and other situation with 
anharmonicity [26]. The Morse Hamiltonian has turned out to have interesting properties in connection with the 
WKB approximation [27l [28| (the WKB quantization gives the exact spectrum), and with supersymetric quantum 
mechanics (29j . 

One can transform the Hamiltonian by means of the displacement operator e lb P so that Ht = e lbp He p — 
^- + A 2 (l — exp[— qi]) 2 . The transformed Hamiltonian is then written in the approximate form (we neglect odd 
powers because of RWA, i.e. we are in the regime A 3> a) 



H = V2\a(N + l/2) + \ 2 (S 1 +S 2 ) (31) 

with 

Following the procedure of the former section we may write 

Sl = e ^L*(-^) ^(N + 1/2) - 1 (33) 

and 

ft- 2 (.*M- i ^)- 5 ^ [ (* + V»-l) (34, 

giving the RWA Hamiltonian 

B .^i! (* + i) + A =„* is (-^) - + tf. (35) 

Using the fact that A ^> a in the RWA validity regime, we expand the Hamiltonian to second order in a/ A 

+ A + (36) 

Note that if we write the (transformed) Morse Hamiltonian as 

H = V2Xa(N + 1/2) + X 2 (l - exph^S]) 2 - (37) 

VV2Aa 2 2V2Aa 



and express the exponentials in a factorized normal form 

exp[a(a + at)] = £ £ " " " ■ (38) 



n=0 k=0 



RWA on the above expression keeps only terms n = k, i.e. the double sum becomes a single sum, leading to Laguerre 
polynomials of order (operator) N. In this form we can recover equation pip via RWA and using equation (137|) . 



V. VALIDITY CHECK OF THE APPROXIMATION 



In this section we numerically study the applicability of the above approximation as funtion of the system param- 
eters. 
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A. Cosine potential 



As pointed out, solutions of the Mathiue equation in closed analytical form do not exist. However, it is well 
known, from Floquet [3(| and Bloch theory [31 j, that the spectrum of a periodic operator in 1-D is determined by 
a continouos quantum number k refered to as quasi momentum and a descrete number n called band index. The 
spectrum is most often represented within the first Brillouin zone [3l| as allowed energy bands seperated by forbidden 
gaps. The characteristics of the spectrum are usually, in a phenomenological way, explained in two different ways; 
starting from a weak potential or a strong one. In the weak limit the particle is moving almost freely in a periodic 
background, which at the degenerate points split the degeneracy. In the oposite, strong potential limit, particles with 
energies smaller than the potential barriers will be quasi bound, and the tunneling rate determines the width of the 
bound energies, band width. As soon as a particle has an energy exceeding the barriers, it will move almost freely. 
Thus, we expect that the ^-dependence is weak only for the lowest quasi bound energy bands when the amplitudes 
go of the cosine potential is large. This is confirmed in fig. [TJ which shows the lowest energy bands of the Mathiue 
equation for q — 1 and g^ = 10. 
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Figure 1: The spectrum E n (k), within the first Brillouin zone, of the Mathiue equation for a potential amplitude g% = 10 and 
wave number 5=1. The cosine function is plotted for clearity, showing that within the "wells" the particle is quasi bound. 

The RWA is only supposed to work in the limit of go 3> q, so it approximate only the lowest quasi bound energies. 
In fig. [2]we present the results of numerical calculations of the error estimete SE(n,go) = \E^ WA — E n \, where E^ w 
is the approximate result for the energy from eq. (p?5)) . 



q 

2ffo 



q 
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1 



(39) 



and E n is the numerically calculated result of the energy by diagobnalization of the truncated Hamiltonian. Here 
the size of the Hamiltonian is 765 x 765 (well within the convergens limits for the eigenvalues). The reason why we 
show the absolute error, and not the relative one, is because for higher values of n, the energies become close to zero 
and the relative error fluctuates greatly in such cases. In the figure, the band index n runs between and 5, hence 
showing the six lowest energies, and it is clear that the approximation breaks down for small couplings go and high 
exitations n as expected. 



B. Morse potential 



As the Morse potential is analytically solvable, no numerical diagonalization of the Hamiltonian is needed. The 
bound energies for the Morse potential are 0, Hfl 



E n = V2Xc 



(n-\ 




a 2 ( 






■\) 










-t r H 





(40) 



Clearly, since A>a, the second anaharmonicity term becomes crucial only for larger excitations n. The anaharmonic- 
ity terms is purely negative resulting in that the highly excited bound states are more densly distributed. Interestingly, 
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Figure 2: The absolute error 5_B(n, g^) = \E^ WA — E n \, between the RWA result and the "exact" numerical result, as a funtion 
of the band index n and the amplitude g$. The approximation is most valid for low excited bands n and strong couplings go- 
Again the wave number q is set to unity. 



from the RWA result of eq. Q36p we note that the second term is positive making the energies to be more sparse for 
high n's. However, this is partly compensated for by the coeffiecient 3/4 infront of the first harmonic part. It should 
be emphasized that only the even terms in the sums (|3l?|) are included since within the RWA odd terms vanishes. 
Thus, it is expected that the method may not be as efficient as for a situation with a purely even potential V(x). It 
is more likely that the obtained eigenvalues approximate the ones for the Hamiltonian with a potential 

f ff(l) -M- 1 + yM!_ 2 yH! (4i) 

Veff[X) A 2 ^ {2k)\ ^ (2k)\ ' [ ' 

Using that e x = ^2 k we expect that the potential (|4Tj) has some kind of "weak" exponential behaviour. In fig. 
[3] examples of the normalized effective potential V e ff(x)/V e ff(l) are given for a — 0.1, 1 and 10. 




-1 -0.5 0.5 



X 

Figure 3: The normalized effective potential V e f f{x)/V e f /(l) defined in eq. (|41|l for a = 0.1, 1 and 10. 



The relative error A E (n,X) = \E^ WA -E n \/E n between the expanded RWA result E^ WA fl3(| and the exact result 
E n (|4"D|) is shown in fig. 2] for the first six eigenvalues and as a function of A. In the plot a = 1, but similar results 
are obtained for other widths a, however, with some increase of A# for larger a. 
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VI. CONCLUSIONS 



We have developed a method to approximate Hamiltonians using a kind of self-RWA. We have shown how some of the 
involved sums may be calculated using symmetrically ordered expressions for Fock-states expectation values of powers 
of the position operator. We have used these expressions to obtain approximations for Hamiltonians corresponding to 
the quantum rotor and the Morse potential. A discussion of the validity of the approximation was given and the 
relation with perturbation theory was explained. The direct link between first order perturbation theory and the 
RWA deepens the understanding of the two methods, and a consequent question would be if higher orders in the RWA 
scheeme could regain higher order perturbation theory. This seems possible, but has turned out to be more subtle 
than expected, mostly because of the non-commutability of the ladder operators. 



Appendix A: CALCULATING SUMS VIA SYMMETRIC ORDER 

In this appendix we use the results of the previous sections and show how one may use it to calculate various sums. 
Like in sectiorllTl no RWA is used and the results are exact. 

We start with some function f(x) which we can expand according to 

k=0 ' fc=0 

Here we have introduced the "symmetric" ladder operators a and a) obeying the regular boson commutator algebra 
[a, aJ] — 1 and which are related to x and p — — ijj as 

£ = 73 ( at + a ) 

(A2) 

P = 75 ^ ~ ' 

Further we have the normalized eigenstates and eigenvalues of the number operator N — a^a; N\n) — n\n), where 
n = 0,l, 2,..., and which in x-basis are 

= (s^) = l^= e -* 2 /iH n {x). (A3) 

As argued in section |ll] the diagonal elements, in the Fock basis, of the function f(x) is identical to the diagonal 
elements of the function f(x,p) where only terms containing an equal number of creation and annihilation operators 
are included. We thus have 

(n\f\n) = (n\f\n). (A4) 
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The L.H.S. can be written, using equation (flO]) . as 



- i ~ /w+")( ) m 

J 2^2 2l l\^ 2J'(2!) a (AT — ?')! 1 j 



and if we express the R.H.S. in the x-basis we obtain 



Below we give some analytically solvable examples. Approximate results is in principle easily achievable, but is left 
out in this paper. 

1. Cosine function 

For a cosine function V(x) = cos(qx) we have V( 2fe 5(0) = (— l) fc <7 2fc and we find [12] in (IA6[) 



(A7) 



t hs -v t- 1 ^ y (-i)V fc _4y (-i)V fc " ! 

• • • Zv 4^1 2 fe (fc!) 2 (n-*)! ° ^ 2 fc (fc!) 2 (n - k)\ 

i <- + OC 2 

R.H.S. = — — -= / cos(gx)e- H 2 (x)dx = e'^ L n (q 2 /2). 



So that 



EwF?^ 1J * (A8) 

Note that closed forms of V(x) = exp(±iqx), V(x) — cosh(qx), V(x) — cosh m (qx) or V(x) = cos m (qx), m — 
0, 1, 2, ... can also be obtained. 

2. Gaussian function 

In the case of V(x) = exp (— a 2 x 2 ), a 2 > 0, the integral in the R.H.S. of eq (|A6|) is analytically solvable [13] 

1 f +0 ° „2= , 2/ ^, „ , . / 2a 2 \" +1/2 a"' ( 2n-l a 2 + 2\ /A . 

e - Q . e -« Hl(x)dx = 2" +1 -^ r — F -n, n; — ; — , (A9) 



2"n!V^7^oo \a 2 + 2j n \ ' ' 2 ' 2a 

where F(.. ...) is the Gauss hypergeometric function. The derivatives are simply y( 2fe '(0) = (— l) k a 2k 
resulting in the L.H.S. 

L ' H ' S ' = h h (2i + 20! 2ifc-)P (n-j)f (A10) 

The above holds also for exponential functions such that > a 2 > — 1, which, when applied to Hamiltonian systems, 
is of interst as they possess an infinite set of bound states. 

3. Hermite polynomials 

When V m (x) — iJ 2 m(^), where m = 0, 1, 2, ... we find [32j 

1 /• + °° 2 2 m / 2 7T7 '?7 1 

R.H.S. = -= / H 2m (x)e- X H 2 (x)dx = 5 — — , (All) 

2»n!V^i-oo 1 ' " V ^ (^!) 2 ( n _ f )!' 1 ' 

while its derivatives are V,n k \o) — 2 ^1)]^' ^2m-2fc(0)- 
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